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Abstract 

The osmotic coefficient of solutions of rod-like polyelectrolytes is con- 
sidered by comparing current theoretical treatments and simulations to 
recent experimental data. The discussion is restricted to the case of mono- 
valent counterions and dilute, salt-free solutions. The classical Poisson- 
Boltzmann solution of the cell model correctly predicts a strong decrease in 
the osmotic coefficient, but upon closer look systematically overestimates 
its value. The contribution of ion-ion-correlations are quantitatively stud- 
ied by MD simulations and the recently proposed DHHC theory. However, 
our comparison with experimental data obtained on synthetic, stiff-chain 
polyelectrolytes shows that correlation effects can only partly explain the 
discrepancy. A quantitative understanding thus requires theoretical ef- 
forts beyond the restricted primitive model of electrolytes. 

1 Introduction 

Polyelectrolytes are macromolecules carrying charged groups which in polar sol- 
vents can dissociate into a highly charged macroion and oppositely charged small 
counterions. 1-5 The high electric field of the macroion strongly couples to the 
counterions, which in turn will tend to partly neutralize the macroion. This phe- 
nomenon has been termed "counterion condensation" . 6 ~ 8 Since polyelectrolytes 
present an ubiquitous class of materials, the quantitative understanding of coun- 
terion condensation is of fundamental importance. 

The strong correlation of the counterions with the macroion reduces their 
thermodynamic activity. Only a certain fraction of ions will be osmotically 
active and contribute to the osmotic pressure n of polyelectrolytes in dilute 
solutions, as Kern was able to show by osmotic measurements as early as 1938. 9 
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The reduced activity can be expressed in terms of the osmotic coefficient 



where Iljd is the ideal osmotic pressure calculated for a solution of counterions 
interacting neither with themselves nor with the macroion. A great number of 
experimental studies of </> for strongly charged polyelectrolytes is now available 
in literature, demonstrating that (p is of the order of 0.2 to 0.3 for univalent 
counterions in the dilute regime. 1-3,10-17 Furthermore, a comparison of exper- 
imental data obtained for different counterions show that specific interactions 
between counterions and the macroion can lead to an additional reduction of 
the measured osmotic pressure. 1 ' 3 In case of monovalent counterions, these spe- 
cific effects are small as compared to the effect of counterion condensation that 
presents the dominant feature of strongly charged polyelectrolytes. A detailed 
discussion of this problem has been given by Mandel. 1 However, no quantitative 
conclusions could be drawn from this comparison so far, since practically all ex- 
perimental data have been obtained from solutions of flexible polyelectrolytes, 
which change their conformation upon lowering the ionic strength. 18 The intri- 
cate connection of the osmotic coefficient with the conformation of the macroion 
is not yet fully understood, therefore studies of flexible polyelectrolytes in solu- 
tion cannot lead to direct conclusions on the validity of theories developed for 
rod-like macroions. 

In this paper we wish to discuss recent theoretical models of rod-like poly- 
electrolytes and compare the calculated osmotic coefficient to results of exper- 
imental investigations of a suitable model system. We restrict the discussion 
to the case of rod-like macroions for which a partial counterion condensation 
occurs, 6 ' 19-21 i.e., part of the counterions remain within a finite distance to 
the macroion and will not be diluted away when lowering the concentration to 
extremely low values. The paper is organized as follows: First the simplified 
situation of the cylindrical cell model and its solution within the framework of 
Poisson-Boltzmann (PB) theory will be discussed. Next we consider theoreti- 
cal refinements designed to overcome limitations inherent to the cell model or 
its mean-field PB solution. In particular, computer simulations will be used 
to assess the importance of correlations neglected in the PB approach. Finally 
recent measurements of the osmotic coefficient for a synthetic rod-like polyelec- 
trolyte will be quantitatively compared to two theoretical models and computer 
simulations. Hence, the present comparison allows to assess to what extend the 
deviations between measurements and PB theory can be accounted for by the 
modern theoretical approaches. 
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2 The cell model 



The cell model is a commonly used way of reducing the complicated many-body 
problem of a polyelectrolyte solution to an effective one-particle theory, i.e., the 
case of a single polyelectrolyte chain in a cell. 3-5 ' 19,20 ' 22 The idea is to partition 
the solution into sub-volumes, each containing only a single macroion together 
with its counterions. Since each sub-volume is electrically neutral, the electric 
field will on average vanish on the cell surface. By virtue of this construction 
different sub-volumes arc electrostatically decoupled to a first approximation. 
One may thus hope to factorize the partition function and reduce the problem 
to the treatment of just one sub-volume, called "cell" . Its shape is assumed to 
reflect the symmetry of the polyelectrolyte. For a solution of rod-like polyelec- 
trolytes with density p p and rod length L this gives a cylindrical cell with the 
radius R being fixed by the condition nR 2 L x p p = 1. The theoretical treat- 
ment is much simpler after neglecting end effects at the cylinder caps. This is 
equivalent to taking L — oo after mapping to the correct density. 

The ions will be described as point-like, but they will electrostatically inter- 
act with the macroion as well as with each other. Their positions may thus be 
strongly correlated. The solvent molecules are not explicitly taken into account; 
rather, they are assumed to form a continuous dielectric background which is 
completely specified by its dielectric constant e r . 

In this simplified situation it can be proved rigorously that the osmotic pres- 
sure is given by the counterion density at the cell boundary times the thermal 
energy fee T. 23 The osmotic coefficient is thus given by the ratio of the ion 
density at the cell boundary divided by the average ion density. 

3 Poisson-Boltzmann theory for the cylindrical 
cell model 

The analytical description taking into account all individual counterions turns 
out to be still too involved, since - as usual - inter-ionic correlations complicate 
matters. The standard way out is neglecting those correlations in a mean-field 
spirit, as is done in the Poisson-Boltzmann description: The ionic degrees of 
freedom are replaced by a cylindrical counterion density, and their interaction 
is approximated by the assumption that the density is locally proportional to 
the Boltzmann- factor, see e.g. Refs.. 19 ~ 21 ' 24 It can be shown that on this mean- 
field level the osmotic coefficient is still given by the ratio between boundary 
density and average density. 25 

Let us briefly recapitulate the Poisson-Boltzmann treatment of the cell model. 
Consider a rod of radius ro and line charge density A > coaxially enclosed in 
a cylinder of radius R. Electroneutrality is achieved by adding the appropriate 
amount of monovalent counterions; no additional salt shall be present. The 
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strength of the electrostatic interactions is conveniently expressed by the Bjer- 
rum length 



,2 
-0 

47T£o£ r fcB7 1 ' 



= (2) 



where eo is the unit charge, e r is the dielectric constant of the solvent, and eq, 
A)b and T have their usual meaning. This definition suggests a dimcnsionless 
way of measuring the line charge density of the rod via the charge parameter £: 

£ = AVe - (3) 

It counts the number of unit charges on the rod per Bjcrrum length and is usually 
called the Manning parameter. In the following only the strongly charged case 
£ > 1 will be discussed. 

Within the framework of PB theory the osmotic coefficient turns out to be 
given by the expression 25 

1 + 7 2 

0PB = ~^P' (4) 

where the dimensionless constant of integration 7 is the solution of the tran- 
scendental equation 

, r 1 i- 1 

7 In — = arctan — h arctan . (5) 

r 7 7 

In the limit of infinite dilution the cell radius R tends to infinity, which implies 
7^0. Note that for 7 — > and £ > 1 the right hand side of Eqn. (5) tends 
to the constant ir, hence the osmotic coefficient as computed by PB theory 
(logarithmically) converges to the well known Manning limiting law (p^, — l/2£. 
At finite densities it is always larger, however. 



4 Beyond the Poisson-Boltzmann solution of the 
cell model 

The theoretical treatment discussed above suffers from various approximations. 
First, the cell model itself is a simplified representation of the polyelectrolyte 
solution: It neglects rod-rod-interactions, it is incapable of describing effects due 
to the finite length of the rods, and it reduces the solvent to a dielectric con- 
tinuum. Second, the mean-field approach discards any inter-ionic correlations 
which can modify the average charge distribution. In this section we briefly out- 
line a few theoretical approaches which try to provide an improved description 
of the physical situation. Our main purpose is to indicate, in which direction 
changes of the osmotic coefficient are to be expected. 
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4.1 Integral equations 

The key approximation in Poisson-Boltzmann theory is the neglect of ionic corre- 
lations by assuming that all iV-particle distribution functions entirely factorize. 
Hence, each ion "sees" only a mean ion density but no individual particles in 
its vicinity. However, the local ionic order around each ion influences its energy. 
Recall for instance that the mean charge density around an ion in a crystal is 
zero, but nevertheless the electrostatic energy per ion is (i) negative, and (ii) 
the Madelung constant depends on the crystal structure. 

Integral equation theories try to explicitly compute higher order correlation 
functions and their impact on the one particle distribution function, i.e., the 
density. The problem of the electric double layer has been treated frequently 
in the past for various geometries. 26-28 What concerns us here is that the 
screening of the macroion by the small countcrions is found to be enhanced. 
Since a larger number of ions close to the rod implies a smaller number of ions 
at the cell boundary, this effect will lower the osmotic coefficient as compared 
to Poisson-Boltzmann theory. 

4.2 Correlation corrected density functionals 

The correlations neglected in the mean field approach can in principle be re- 
captured by adding a correction to the free energy functional corresponding 
to Poisson-Boltzmann theory. 29-31 A local and stable functional based essen- 
tially on a Debye-Hiickel solution of the one component plasma has recently 
been suggested by Barbosa et at. 32 It leads to an enhanced counterion con- 
densation in the vicinity of the macroion in very good quantitative agreement 
with Molecular Dynamics simulations. This again implies a lower density at 
the outer cell radius and thus a smaller osmotic coefficient. The extended free 
energy functional translates to a differential equation more complicated than 
the PB equation. However, the equilibrium ion profiles can be conveniently de- 
termined by directly implementing the functional minimization problem along 
the lines suggested in Ref.. 33 A comparison of Poisson-Boltzmann theory, com- 
puter simulations and experiments with this correlation corrected theory will 
be presented below. 

4.3 Debye-Hiickel-Bjerrum theory 

An attempt to go beyond the cell model with the aim of incorporating rod-rod- 
interactions has been suggested by Kuhn et al.. 34 A solution of charged rods is 
investigated within the Debye-Huckel-Bjerrum theory, which uses the linearized 
form of the Poisson-Boltzmann equation to solve for the charge distributions. 
The neglected nonlinearities are approximately recaptured by the explicit intro- 
duction of clusters consisting of a rod and condensed countcrions. Minimizing 
the free energy of the system, which is constructed as a sum of electrostatic 
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and entropic terms, gives the cluster distribution and thermodynamic quanti- 
ties like the osmotic pressure. No cell-description is used, rather, the interaction 
between individual rods is explicitly computed and their mutual orientations is 
averaged over. For the computation of the screened rod-potential (and at a few 
other stages) end effects are neglected, which in the dilute limit entails the Man- 
ning limiting behavior of infinite rods. For the case of added salt the theory 
predicts an osmotic coefficient even below the Manning predictions, but it is 
unfortunately not straightforward to extend it to the salt-free case. A simpler 
version for the salt free case 35 that neglects rod-rod interactions, again shows 
the osmotic coefficient to be always below the Manning limiting value. This, and 
the fact that it decreases upon increasing polyelectrolyte concentration, raises 
concerns as to whether the employed approximations lead to a well-controlled 
improvement beyond the cell model. 

4.4 Finite rods 

An alternative theory beyond the cell model has been put forward by Nyquist 
et al.. m Here the aim is to explicitly account for the finite length of the rods. 
It is argued that in the limit of infinite dilution a finite rod will not produce 
Manning condensation and the ions will become free. This implies the correct 
limiting osmotic coefficient of 1 instead of l/2£ at zero density. At finite density 
the condensation is approximated by a two state model (ions are either free 
or condensed, and each rod carries the same fraction of condensed ions). The 
partition function is computed in a random phase approximation, which finally 
gives the equilibrium condensed fraction as a function of density. In the dilute 
regime this fraction is found to increase with density and the osmotic coefficient 
is correspondingly found to decrease. Furthermore, the calculated osmotic co- 
efficient is larger than the Poisson-Boltzmann result. This is in strong contrast 
to any solution of the problem which assumes infinite rods. 

Up to now there is no systematic comparison of theory and experiment, 
however, the model seems to overestimate considerably the osmotic coefficient 
when compared to experimental data. 10,17 Possible reasons for this shortcoming 
may be related to the two state model or the random phase approximation 
employed in Ref.. 36 

5 Simulations 

After restricting to the simplified models mentioned above, simulations can be 
very successfully used to test the available theories, demarcate their range of 
applicability, and indicate the trends of deviations which are to be expected. 
Consequently, versions of the cylindrical cell model have been used several times 
in the past as a basis for simulations. 24,37-39 An important observation has been 
that condensation is generally enhanced compared to Poisson-Boltzmann theory, 



6 



which is frequently attributed to the correlations neglected by the latter. 

In this article we use molecular dynamics simulations to determine the os- 
motic coefficient for a cell model, the parameters of which have been mapped 
to a specific polyelectrolytc: poly(para-phenylene) (PPP, see Fig. 1 and below). 
We take a cubic simulation box of length L^, one charged rod parallel to an 
edge, and the correct number of monovalent counterions which leave the system 
electrostatically neutral. The ions interact via a purely repulsive Lennard-Jones 
potential, 18 giving them a diameter of 4.4 A. For the rod-ion-interaction a sim- 
ilar potential was used, in which the hard core was shifted to larger radii, such 
that the distance of closest approach between an ion and the rod was 7 A. This 
value has been inferred from an analysis of small angle (neutron or X-ray) scat- 
tering experiments on ppp. 4 °- 42 Each ion carries a negative unit charge in its 
center and the rod carries a sequence of positive unit charges along its axis at a 
distance of 2.15 A. Together with a Bjerrum length of 7.31 A (corresponding to 
water at 40 C) this gives a charge parameter £ = 3.4. After switching on peri- 
odic boundary conditions this geometry yields an infinite array of infinitely long 
rods sitting on a square lattice. When comparing to the cylindrical cell model 
we will map the rectangular cell belonging to one such rod with a cylindrical 
cell of radius R = Lb/y/n. 43 

The electrostatic interactions in this periodic boundary geometry were com- 
puted with the help of P 3 M routines, 44 and a Langevin thermostat 45 combined 
with a velocity- Verlet-integrator 46 (with time-step 0.01 in Lennard-Jones units) 
was implemented to drive the system into the canonical state. The saturation of 
the electrostatic energy was used to test for equilibration. A more detailed de- 
scription can be found in. 47 We simulated the rods at a monovalent counterion 
concentration of 2.18, 6.79, and 10.2 mmol/1, with a number of 4.5, 1.8, and 1.2 
million MD steps, respectively. The corresponding number of monovalent ions 
in the simulation cell were 277, 1256, and 1024. The osmotic pressure was taken 
as the average of the xx- and yy-component of the stress tensor, given that the 
rod points along the z-axis. 

6 Experiment 

From its derivation it is obvious that the cylindrical cell-model can be valid 
only if the macroions are sufficiently stretched. DNA in salt-free solution can 
be regarded as such a model polyelectrolyte. 10,41 ' 48 However, the stability of 
the helical conformation is severely impeded if the ionic strength is low, and 
measurements in salt-free solution must be carried out at low temperature to 
avoid the melting of the helix. To the author's best knowledge the study of Auer 
and Alexandrowicz 10 has been the only one in which strictly salt-free solutions 
of DNA have been studied. A more recent investigation employed solutions 
of DNA with additional 2mM and 10 mM added salt. 48 However, in this case 
counterions and salt ions will establish a Donnan equilibrium, 2 ' 3 ' 49 which results 
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in a drastic reduction of the small-ion contribution to the pressure. Its dominant 
origin at low polyelectrolyte density is not yet fully understood. 50 

Recently, an investigation of the osmotic coefficient of a synthetic rod-like 
polyelectrolyte has been presented. 17 The system studied there consists of a 
poly(p-phenylene) (PPP) backbone to which positively charged side-chains have 
been attached (see Fig. 1). The fully aromatic backbone exhibits a high stiff- 
ness with a persistence length of approximately 22 nm. 51 The number-average 
contour length of PPP-macroions used in the experiment is of the same order. 
Hence, the PPP-polyelectrolytes may be treated as rod-like in very good ap- 
proximation. The system has a charge parameter £ = 3.4 of similar magnitude 
as DNA (£ = 4.2). The excellent chemical stability of the PPP-chains allows to 
study these polyelectrolytcs in salt free solutions without any problem. 

In Ref. 17 the osmotic coefficient of the PPP system was measured for two 
kind of counterions, iodine and chlorine. Both sets of experimental data, which 
have been obtained in virtually salt free solution, clearly indicate that the 
Poisson-Boltzmann solution of the cell-model makes a fairly good prediction 
for the osmotic coefficient, but systematically overestimates it upon closer look. 
This is also in agreement with the earlier results obtained on DNA by Auer and 
Alexandrowicz. 10 Moreover, it became apparent in this study that there are spe- 
cific interactions of the counterions with the macroion, which can be as large as 
the differences from the Poisson-Boltzmann solution itself. This was concluded 
from the fact that chlorine counterions lead to a considerably smaller (f> than 
iodine counterions. Such specific interactions have been found for many flexible 
polyelectrolytes and extensively discussed in the earlier literature. 1 ' 3 The errors 
in the experiment are hard to estimate, but are on the order of 7%. 

7 Comparison of theory, simulations and exper- 
iment 

It has to be understood that an analysis of the osmotic coefficient </> measured in 
salt free solutions within the framework of the Poisson-Boltzmann solution of the 
cylindrical cell model proceeds virtually without adjustable parameters. Here, 
(f> is solely determined by the charge parameter £, which is fixed by chemistry, 
the rod radius ro, which has been measured experimentally, and the polyelec- 
trolyte concentration, which is fixed by the experiment. The latter parameter 
determines the cell radius R (see the discussion in Ref. 17 ). For the theoretical 
or simulational modeling those values were taken as input parameters and were 
not adjusted such as to fit the measured data. 

Fig. 9 summarizes the results. It shows the osmotic coefficient <j> of a PPP 
solution as a function of counterion concentration c c as predicted by Poisson- 
Boltzmann theory, the DHHC correlation-corrected treatment from Sec. 4.2, 
Molecular Dynamics simulations and experiment. 
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To begin with, it is important to note the fine vertical scale. Poisson- 
Boltzmann theory predicts <f) to be smaller than 1 and vary roughly within 
the range 0.18 . . . 0.22. The measured values accumulate around 0.18 (iodine) 
and 0.16 (chlorine). Hence, the dominant change in <fr, a reduction by a factor 
of 5, is correctly accounted for. However, on the enlarged scale of Fig. 9 it is 
visible that the measured values are systematically lower than the prediction, 
although still higher than the Manning limit l/2£ of infinite dilution. 

Both the correlation-corrected DHHC theory as well as the simulations, 
which in principle capture all kinds of ion correlations, show a decrease in the 
osmotic coefficient. Since these two totally different approaches agree so well, 
we believe that they indeed give a good description of the influence of correla- 
tions. However, they do not lower the osmotic coefficient sufficiently enough as 
to agree with the experimental results. 

The deviation from the Poisson-Boltzmann curve increases for higher den- 
sities, which is true for the DHHC and simulational part as well as for the 
experiment. This appears plausible if one recalls that correlations become more 
important at higher densities. Unfortunately, the increasing systematic experi- 
mental errors in this regime are very difficult to quantify. Due to this fact, we 
give no significance to the apparent slight decrease of <fi for the I - ions. 

8 Discussion 

The fact that Poisson-Boltzmann theory overestimates the osmotic coefficient 
as been observed previously. Careful studies of typical flexible polyelectrolytes 
in solution (see Ref. 1 ' 3 and further references given there) indicated that agree- 
ment of the Poisson-Boltzmann cell model and experimental data could only be 
approached if the charge parameter £ was rcnormalizcd to a higher value. The 
motivation given for this ad hoc modification was the assumption of a locally 
helical or wiggly main chain. Hence, the counterions "see" more charges per 
unit length, i.e., a macroion having a higher charge parameter. However, the 
results obtained for stiff-chain macroions 17 show that the osmotic coefficient 
is lower than the Poisson-Boltzmann results even for systems where the local 
conformation of the macroion is absolutely rod-like, therefore this explanation 
can not be used. Also, the measured as a function of density has a functional 
form different from the Poisson-Boltzmann prediction, and no value of £ gives 
a curve that fits all data points. 

Since Poisson-Boltzmann theory neglects all ion-ion correlations, it is tempt- 
ing to assume that their incorporation into the theoretical treatment would re- 
solve the discrepancy. However, our comparison shows that these correlational 
effects can be made responsible only for part of the deviations. Since the dif- 
ferent approaches using a correlation-corrected density functional theory and 
Molecular Dynamics simulations agree very well with each other, we believe 
that the discrepancy between them and the experiment is not due to theoretical 
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or simulational errors, but has to be looked for in the approximations underlying 
the cell model itself. 

If there were any excess salt present in our system, the osmotic coefficient 
would be reduced. 3,48 ' 52 However, we performed numerical solutions of the 
Poisson-Boltzmann equation, taking into account various amounts of excess salt 
as well as the constraints due to the Donnan equilibrium, which showed two 
things (see Fig. 9): First, the functional form of <j>(c c ) is not compatible with a 
roughly constant value of the osmotic coefficient. Rather, it approaches the salt- 
free result at high c c and drops strongly (below the Manning limit) at low c c . 
Second, our measurements at small counterion concentration are incompatible 
with the assumption of an excess salt concentration exceeding a few micromolar. 
In fact, the strong decrease of the osmotic coefficient at low concentrations, 
even below the Manning limiting value for the Cl _ ions, can be explained by 
different very small excess salt concentrations for both experimental sequences 
(2 mikromol for I , 20 mikromol for CI - , compare Fig. 2). However, at our 
high concentrations, where the discrepancy is most apparent, any remaining 
excess salt is least relevant, and hence cannot explain neither the difference in 
the osmotic coefficient between Cl~ and I - ions nor the decrease of <p compared 
to the PB prediction. 

One obvious approximation of the cell model is the neglect of end effects 
of the rods by assuming them to be infinitely long. However, as we have dis- 
cussed in Sec. 4.4, correcting this point would lead to an increased osmotic 
coefficient, since finite rods show less condensation (indeed, none in the limit of 
zero density). Moreover, the theoretical treatment of finite rods by Nyquist et 
al. 36 produces too large coefficients even at finite densities. Hence, finite size 
effects of the rods are unable to explain the additional reduction of the osmotic 
coefficient. 

The theoretical approach mentioned in Sec. 4.3 yields a reduced osmotic co- 
efficient. However, it is unable to predict values above the Manning limit, and 
shows a decrease of </> with increasing concentration, contrary to the experimen- 
tal indications. The slight apparent decrease for the iodine is, if the errors are 
considered, not significant. 

The remaining deficiency of the cell model is the neglect of any molecular de- 
tail, in particular: The solvent is described by a dielectric continuum. However, 
hydration effects are well known to be important in many circumstances 53-55 
and to depend on the particular objects to be hydrated. 55 This of course applies 
to the counterions as well as to the charged groups of the macroion. Indeed, the 
difference between the osmotic coefficient measured for chlorine and iodine ions 
indicates that they are not just charged spheres. Similarly, the macroion is not 
just a charged rod. 

Unfortunately, those last problems with the cell model are the hardest to 
deal with theoretically. Even predicting the direction into which changes are 
to be expected are very complicated. 56 Obviously, much work still has to be 
done in order to arrive at an improved quantitative prediction of the osmotic 
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coefficient of rodlike polyelectrolytes, since all the "obvious" modifications will 
not suffice, as we tried to show. 

9 Conclusion 

We tested the ability of the cell-model to explain the measured osmotic co- 
efficient of solutions containing rodlike polyelectrolytes and monovalent coun- 
terions. We compared the Poisson-Boltzmann solution of the cell model, an 
improved local density approximation (DHHC), and simulational results of the 
model to recent experimental data. 

Our findings are that the Poisson-Boltzmann solution of the cell model sys- 
tematically overestimates the osmotic pressure. We showed that only a part of 
this discrepancy between theory and experiment is due to the neglect of cor- 
relations in the mean field approach, which will lower the osmotic pressure. 
We argued that other simple explanations like a locally curved main-chain con- 
formation or the neglect of remaining excess salt do not apply here. We also 
outlined why dropping the assumption of infinite rods would in fact worsen the 
disagreement 

We are led to the conclusion that the remaining discrepancy between exper- 
iment and a theoretical description of the system should be looked for on the 
level of molecular detail. This includes a better description of the solvent, hy- 
dration effects, or van der Waals forces. We are well aware of the fact that this 
requires formidable theoretical or numerical efforts, but we would like to stress 
that a full agreement between theory and experiment cannot and should not be 
expected on the level of a restricted primitive model for electrolytes, even if all 
correlations are properly taken into account. We believe that our comparison 
provides some evidence to this statement and thereby hope to motivate future 
theoretical and experimental work to resolve these issues. 
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Figure 1: Constitution formula for poly(p-phenylene). The fully aromatic back- 
bone exhibits an excellent chemical stability. Its uncharged precursor has a 
persistence length of approximately 22 nm. The degree of polymerization used 
in the studies in Ref. 17 ' 42 was located between 20 and 40. Therefore the contour 
length equals approximately one persistence length at most. 
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Figure 2: Osmotic coefficient as a function of countcrion concentration c c for 
the poly(p-phenylene) systems described in the text. The solid curve is the PB 
prediction of the cylindrical cell-model, the dashed curve is due to the correlation 
corrected PB theory from Rcf. 32 (see Sec. 4.2). The circles are experimental 
results using iodine (full) and chlorine (empty) counterions. Errors are of the 
order of 7% for all values, but only indicated for one measurement. The crosses 
originate from the MD simulations described in the text. Furthermore, the 
thin dashed lines are PB predictions assuming excess salt concentrations of 2 
micromolar (upper) and 20 micromolar (lower), which are to be understood as 
crude estimates of the actual excess salt content. The Manning limiting value 
of l/2£ is also indicated as a horizontal dashed line. 
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